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We provide a theory of the deflection of polar and non-polar rotating molecules by inhomogeneous 
static electric field. Rainbow-like features in the angular distribution of the scattered molecules are 
analyzed in detail. Furthermore, we demonstrate that one may efficiently control the deflection pro- 
cess with the help of short and strong femtosecond laser pulses. In particular the deflection process 
may by turned-off by a proper excitation, and the angular dispersion of the deflected molecules can 
be substantially reduced. We study the problem both classically and quantum mechanically, taking 
into account the effects of strong deflecting fleld on the molecular rotations. In both treatments we 
arrive at the same conclusions. The suggested control scheme paves the way for many applications 
involving molecular focusing, guiding, and trapping by inhomogeneous flelds. 
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I. INTRODUCTION 



Deflection of molecules by inhomogeneous external 
fields is an important subject of molecular physics, which 
continues to attract a lot of attention in recent years 
1-7|. The external fields can be magnetic or electric 
1-4 1, or even optical fields of strong lasers [sl-fl^. By 
controlling molecular translational motion with external 
fields, novel elements of molecular optics can be real- 
ized, including molecular lens 0, Q and molecular prism 
0. Defiection by external fields is also used as a tool to 
measure molecular polarizability |3] and molecular dipole 
moment. The mechanism of molecular deflection by a 
nonuniform static electric field is rather clear. For a 
non-polar molecule, the field induces molecular polariza- 
tion, interacts with it, and deflects the molecules along 
the interaction energy gradient. For a polar molecule, 
the field interacts with the molecular permanent dipole 
moment as well. As most molecules have anisotropic 
polarizability or/and a permanent dipole moment, the 
deflecting force depends on the molecular orientation 
with respect to the deflecting field. Previous studies 
on molecular defiection have mostly considered randomly 
oriented molecules, for which the defiection angle is some- 
how dispersed around the mean value determined by the 
orientation-averaged polarizability (or dipole moment). 
This dispersion was observed via broadening of the scat- 
tered molecular beam, and was re por ted in the "two- 
wire" electric field experiments [3. Il3l [l3 and also in 
the multipole electric field experiments [15l - [l7| . More re- 
cently, this kind of rotation-induced dispersion in molec- 
ular scattering by static electric fields was used as a se- 
lection tool in experiments on laser-induced molecular 
alignment 18|. The field- molecule interactions become 
intensity-dependent for strong enough fields due to the 
field-induced modification of the molecular angular mo- 
tion [l^, [iO] • This adds a new ingredient for controlling 
molecular trajectories d, [l^, Ho, HH . 

Recently, we showed that molecular deflection by 
strong fields of focused laser beams can be significantly 
affected and controlled by pre-shaping molecular angu- 



lar distribution before the molecules enter the interac- 
tion zone [12]. This can be done with the help of nu- 
merous recent techniques for laser molecular alignment, 
which use single or multiple short laser pulses (transform- 
limited, or shaped) to align molecular axes along certain 
directions. Short laser pulses excite rotational wavepack- 
ets, which results in a considerable transient molecular 
alignment after the laser pulse is over, i.e. at field-free 
conditions (for reviews on field-free alignment, see, e.g. 
[13, [13) • Field- free alignment was observed both for 
small diatomic molecules as well as for more complex 
molecules, for which full three-dimensional control was 
realized [gil^. 

In the present paper we extend this approach to molec- 
ular deflection in static electric fields, and demonstrate 
that the average scattering angle of deflected molecules 
and its distribution may be dramatically modified by a 
proper field-free pre-alignment with ultra-fast lasers. An 
important difference between the scattering in a static 
field and an optical field is due to the role of the 
molecular permanent dipole moment. Dipole interac- 
tion with the laser beams averages to zero because of 
the fast oscillations of the optical fields, however this 
kind of interaction becomes dominant for polar molecules 
placed in the static electric fields. In the present pa- 
per, we analyze in detail interaction of rotating molecules 
with the static electric fields (taking into account both 
the dipole- type and polarization- type interactions), and 
demonstrate that laser-induced pre-alignment provides 
a fiexible tool for controlling molecular motion in these 
fields. 

In Sec. [TT]we present the deflection scheme, and pro- 
vide heuristic arguments on the anticipated role of molec- 
ular rotation on the scattering process (both for thermal 
molecules and molecules pre- aligned by additional laser 
pulses). In Sec. IIIII we provide a full classical treatment 
of the problem, which supports the heuristic predictions 
of Sec. ini An alternative classical analysis based on the 
formalism of adiabatic invariants is given in Sec. I VII and 
it leads to the same results. In Sec. |lV]we support our 
control approach by means of a full-scale quantum me- 
chanical analysis. Finally, we summarize our results in 
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Sec. m 



II. MOLECULAR DEFLECTION 

We consider a deflection scheme that is based on the 
interaction between a hnear molecule and an inhomo- 
geneous static electric field. In particular, we follow 
the lines of experiments [1, [2^, in which a collimated 
particle beam goes through a long deflector made of two 
cylindrical electrodes' faces. Electrical field F between 
the t wo p oles is equivalent to an "electrical two- wire 
field" [29|. This geometry allows to obtain an electric 
field F and a field gradient dF/dz which are nearly 
constant over the width of the collimated molecular 
beam (see Fig. [T]). 




V, 
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FIG. 1: The deflection scheme. Linear molecules, initially 
moving in the x direction (with velocity Vx), enter a static 
electric fleld (directed along the z axis). They are deflected 
by the fleld gradient, and get the deflection velocity v^. 

The interaction potential of a linear molecule in the 
static field is given by: 



U = --F^ {Aacos^9 + a±) - nFcosB, (1) 

where F is the electric field; ay and a± are the compo- 
nents of the molecular polarizability along the molecular 
axis, and perpendicular to it, respectively, and fj, is the 
permanent dipole moment. Here 6 is the angle between 
the electric field direction (along the laboratory z axis) 
and the molecular axis. A molecule initially moving along 
the X direction acquires a velocity component Vz along 
z-direction. We consider the perturbation regime cor- 
responding to a small deflection angle, 7 w v^/vx and, 
therefore, assume the molecules are subject to the flxed 
values of the field and field gradient {F and Vi^, respec- 
tively) inside the defiector. 

The deflection velocity is given by: 



(2) 



where M is the mass of the molecules. The time- 
dependence of the force (and of the potential U) in 



Eq.([2]) comes from two sources: projectile motion of the 
molecule through the deflector, and time variation of the 
angle 9 due to molecular rotation. For simplicity, we 
neglect the edge effects at the entrance and exit of the 
deflector. 

Since the rotational time scale is the shortest one in 
the problem, we average the force over the fast rotation, 
and arrive at the following expression for the deflection 
angle, 7 = Vx/v^: 



{F [a\\A2 



(I-A2 



Mvx ' 



(3) 



Here Ai^2 = cos^''^9 denotes the time-averaged value of 
cos^^^ 9, and td is the passage time through the deflector. 
The quantities Ai^2 depend on the relative orientation of 
the vector of angular momentum and the direction of the 
deflecting field. It is different for different molecules of 
the incident ensemble, which leads to the randomization 
of the deflection process. 

We provide below some heuristic classical arguments 
on the anticipated statistical properties of Ai^2 in the 
case of weak fields that do not disturb significantly the 
molecular rotation. We start with the simplest case of 
a linear molecule with fi — 0, which rotates freely in 
a plane perpendicular to the vector ~f of the angular 
momentum (see Fig. ([2])). 




FIG. 2: A molecule rotates with a given angular momentum 
J that is randomly oriented in space. 9j is the angle between 
the angular momentum and the laboratory z axis. 

The projection of the molecular axis on the vertical 
z-direction is given by: 



cos 9{t) — cos{ujt)sm9j, 



where 9j is the angle between J and z-axis 
angular frequency of molecular rotation. 
By averaging Eq. ^ over time, we obtain 



(4) 



and w is the 



Ai = cos 6* = 



A2 = cos2 ^' = 2 



(5) 



For random and isotropic orientation of vector J 
in space, the probability density for 9j distribution is 
l/2sin(6'j). The mean value of the deflection angle is 
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then (7) = 7o, where the constant 70 presents the aver- 
age deflection angle for an isotropic molecular ensemble: 



7o 



FWFtg 

MVr 



(6) 



Eq.® allows us to obtain the distribution function, 
f{A2) for A2 (and the related deflection angle) from the 
known isotropic distribution for 9j. Since the inverse 
function 0j{A2) is multivalued, one obtains 



2 1 



i=l 



1 



(7) 



where we summed over the two branches of 9j{A2)- This 
formula predicts an unimodal rainbow singularity in the 
distribution of the scattering angles at the maximal value 
7 = 7o(q;|| +a±)/2a (for A2 = 1/2), and a flat step near 
the minimal one 7 — "fQa±/a (for A2 — 0). These 
results are similar to the ones derived by us previously 
for molecular scattering by oscillating optical fields [23 |. 

As the next example, we consider the opposite case of 
polar molecules with /i 7^ and negligible polarization- 
type interaction. For the sake of simplicity, we start with 
a 2D model, i.e. for a molecule that rotates with no 
azimuthal momentum. In the limit of E/fiF <^ I (E 
is the rotational energy), the molecular axis is trapped 
by the electric field, and Ai ^ 1. As £' is increased, 
the molecules may still be trapped, but Ai < 1 and for 
trapped molecules with high enough E, A\ < 0. The lat- 
ter happens because the molecules spend most of the time 
being against the electric field when performing nonlin- 
ear angular oscillations. As the energy is increased even 
more, the molecules become untrapped, and perform full 
rotations. In this case, we expect Ai < due to the same 
reason: the molecules accelerate their rotation when the 
dipole moment tends to be parallel to the electric field, 
and they decelerate it when the dipole moment looks 
against the field. As a result, the time-averaged value 
of cos 9 is negative. Considering a 2D molecular rotation 
in the presence of the electric field, we can write: 



dt 



dO 



2 y/E + JlFcost 



(8) 



where / is the moment of inertia. Assuming the un- 
tra ppe d regime {fj,F/E < 1), the rotation period is given 

by 



T, 



period 



do 



81- 



y/E + /iF COS ( 

2 ,n , 



Ve + JIf '2 

The time averaged value of cos 9 is 



(9) 



Ai 



period 

E 



d9 



y/E + fiF cos 9 



cos 9 



(10) 



Here r 



E+fj,F 



E (^^,r) and (f , J') are the first and 
second order elliptic integrals, respectively. From Eq. 
([TI?|) we learn that Ai ~ —fiF/AE in the limit of weak 
fields, fJ.F/E ^ 1. This value is negative, as discussed 
above. 

Summarizing, in the 2D approximation, Ai = 1 when 
E/fiF < 1, it is close to when fiF/E < 1, and it takes 
negative values in-between. This assumes the existence 
of a negative minimum of Ai as a function of E / jj,F . 

The properties of the A2 quantity are somehow differ- 
ent in the considered 2D model. For low E, the molecular 
rotation is suppressed, and A2 = 1. The period of angu- 
lar oscillations of the trapped molecules {fiF/E > 1) is 
given by: 



T, 



period 




d9 



y/E + fiF cos 9 



'-f T- - 

2\j fiF \2' r 
The time-averaged value of cos^ 9 is: 



(11) 



A2 = 



Tperiod 

2 E 



"'(-tS^) d9cos^9 



y/E + iiF cos 9 



4 E g(f,^) 
S^^FFU,^] 



(12) 



where r, £^ (f , 7 



and F (f , i) were defined above. 
The function in Eq. (|12p has a local minimum at A2 = 
0.279, which suggests a rainbow peak in the distribution 
of A2 in the case of a smooth distribution of the parame- 
ter E/fj,F. For high enough energy, the molecules rotate 
almost as free rotors, and we expect a rainbow peak in 
the A2 distribution A2 — 0.5, as was suggested by 
Eq.©. 

Finally, in order to complete our analysis, we consider 
the field-affected molecular rotation in 3D case, find nu- 
merically and plot the quantities of Ai^2 (FigS- El and 
m respectively) for different values of dimensionless ro- 
tational energy and azimuthal canonical momentum (the 
details of the calculations can be found in the next sec- 
tion). For the lowest possible negative values of the to- 
tal energy, E both plots demonstrate angular trapping 
(^1 — 1,A2 ~ 1). For small values of the azimuthal mo- 
mentum we observe a negative shift (around E « fj,F) 
of the maximum of Ai, the nature of which has been 
already discussed. In this limit, we also observe a mini- 
mum at A2 = 0.274 and a peak at A2 = 0.5, which is in 
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agreement with the previous 2D model. For high rota- 
tional energies (and high azimuthal momentum), Ai ~ 0. 
Finally, for high azimuthal momentum we observe also 
a strong decrease in the A2 values, since the molecules 
mainly rotate in the xy plane in this limit. 



0.5 




FIG. 3: Contour plot of ^1 for different values of dimen- 
sionless rotational energy (i.e. h = vertical axis) and 

azimuthal canonical momentum (i.e. I2 = horizon- 
tal axis). The blank part of the figure corresponds to non- 
physical combinations of the parameter values. 




FIG. 4: Contour plot of ^2 for different values of h and I2 {li,2 
were defined in the caption of Fig. [3|. As in Fig. |3] the blank 
part of the figure corresponds to non-physical combinations 
of the parameter values. For low azimuthal momentum, a 
minimum at A2 = 0.274 and a peak at A2 = 0.5 are observed, 
which suggests two rainbows in the A2 distribution. 



III. CLASSICAL TREATMENT 

Consider a classical rigid rotor (linear molecule) de- 
scribed by the Lagrangian: 



P sin^ d 



2 



^2 ( Aa cos^ e 
fiF cos 9, 



(13) 



where 6 and 4> are Euler angles, and / is the moment of 
inertia. The canonical momentum for the (j) angle 

P^^I^sin^e (14) 

is a constant of motion as is a cyclic coordinate. The 
canonical momentum Pg is given by 

Pe = 10. (15) 
The Euler-Lagrange equation for the 9 variable is 
d dL dL 



which leads to 
d^9 ^ " ^ ^ 



dt^ 



Pa 



cos 9 F'^Ac 



■ sm 



sm 



fiF 

I cos f — — sm 



^sin^. (17) 



When considering a thermal ensemble of molecules, it 
is convenient to switch to dimensionless variables, in 
which the canonical momenta are measured in the units 
of pth = I^th^ with LUtfi = y^ksT/I, where T is the tem- 
perature [3l|, and is the Boltzmann's constant. By 
setting = P^/pth, P'g = Pe/Pth. and t' = ujtht, Eq.((I7D 
becomes: 

"^'^ " ■ " " Dsm9 (18) 



dt'^ 



PL 



sm^9 



C sin 9 cos 9 



where C = F'^Aa/ikeT), and D = fiF/{kBT). 

Considering a deflecting field that is adiabatically in- 
creasing to its final value F (adiabatic with respect to 
the molecular rotational dynamics) , we numerically solve 
Eq. ([T5)) and find the time dependent values of cos d{t) and 
cos"^ 9{t). 

In order to find the Ai^2, we calculate: 



A,2(i) 



1 



s^'^ 9dL 



(19) 



and consider Ai^2 = cos^'^ 9 as the limit value to which 
Ai,2{t) converges as t — tp 00. Here tp is the rising 
time in which the defiecting field reaches its maximal 
value F. 

The probability distribution of Ai,2 is given by: 



/(A,2) 



d9{0)d(j){0)dPl){0)dP;i,{0) 



X S{Ai,2 - cos^ '^ 9) 

X /(^^(O),0(O),P^(O),F;(O)), 



(20) 
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where 



, 1 
/ = -5-^ exp 




(21) 



is the thermal distribution function. 

Deflection of thermal molecules. Using this ap- 
proach, we considered distribution functions for Ai^2 
(and corresponding distributions of the deflection angle 
7) for a thermal beam of KCl molecules. For the chosen 
rotational temperature T = A.GSK, the typical "ther- 
mal" value of the angular momentum is Jt = 5, where 
Jt = ^JksT /{hBrc), Br is the rotational constant and c 
is the speed of light. We plot the distribution functions 
for moderate {1.8-lO^V/m) and strong {1.8-lO'^V/m) val- 
ues of the deflecting field at Figs. [S] and [51 respectively. 
For the moderate field, C = 5.81 • 10"^ and D = 0.96, so 
that the dipole-field interaction is comparable with the 
typical thermal rotational energy, while the polarization- 
type interaction is negligible. In this case, a sizable por- 
tion of molecules are trapped by the field, which is re- 
flected in the high positive values of Ai (Fig. [S^). The 
untrapped molecules are performing full rotations, and 
they contribute to the negative values of Ai- Fig. [5)d 
presents the distribution of A2, and it shows two rain- 
bows (the first one at approximately 0.28, and the sec- 
ond one at 0.5), as is expected from the discussion in the 
previous section. The distribution of the deflection an- 
gles (not shown here) is similar to the distribution of Ai , 
as the contribution from the polarization-type interac- 
tion (proportional to A2) is negligible in the considered 
numerical example. The A2 distribution for the KCl 
molecule may be directly measured in a deflection exper- 
iment that combines a homogeneous static field and an 
inhomogeneous laser field. The static field will define the 
distributins of Ai and A2 , and the laser field will defiect 
the molecules according to the A2 values j23 |. 

In the case of a strong field (1.8 • lO'^V/m) shown at 
Fig. m £) = 9.62, and C is stfll negligible. Now the 
dominant portion of molecules is highly trapped by the 
electric field, and the distribution of ^1 is shifted accord- 
ingly to the higher positive values (Fig. [6^). The rainbow 
at 0.5 in the A2 distribution (Fig. [SJj) practically disap- 
pears due to the increased amount of molecules with the 
suppressed rotation. 

Finally, in Fig. [7] we present the case of a strong de- 
fiecting field, but at higher temperature. Though the 
field is as strong as in Fig. [Bl D value is similar to the 
one of Fig. [5l and the curves are correspondingly similar. 

Deflection of pre-aligned molecules. Assume now 
that the molecules are subject to a femtosecond pre- 
aligning pulse polarized in z-direction at i = 0, before 
they enter the defiecting field. The interaction of molec- 
ular permanent dipole moment with the laser pulse aver- 
ages to zero because of the fast optical oscillations. The 
polarization-type interaction is given by the first term 
in Eq. (IT|), in which F is replaced by the envelope e 
of the femtosecond pulse, and an additional factor of 
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FIG. 5: Statistical distributions of (a) Ai and (b) A2 for a 
thermal beam of KCl molecules in a moderate deflecting field 
(1.8 ■ lOV/m). Jt = 5; C = 5.81 • 10"'* ; D = 0.96 
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FIG. 6: Statistical distributions of (a) Ai and (b) A2 for a 
thermal beam of KCl molecules in a strong deflecting field 
(1.8 ■ WW/m). Jt = 5; C = 5.81 • 10"'' ; D = 9.62. 



0.5 is added due to the oscillatory nature of the opti- 
cal field. We assume that the pulse is short compared 
to the rotational period of the molecules, and consider it 
as a delta-pulse. The rotational dynamics of the laser- 
kicked molecules is then described by the same formalism 
as above, but with Pg{0) replaced by 



P^(O)^P^(O)-F,'sin(20(O)). 



(22) 



Here P^. = Ph/ ^ksTI is a properly normalized kick 
strength of the laser pulse, with P given by: 



P = {I /Ah) •(«,,- a^) 



\t)dt. 



(23) 



Here we assumed the vertical polarization (along z-axis) 
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FIG. 7: Statistical distributions of (a) Ai and (b) A2 for a 
thermal beam of KCl molecules at higher temperature that 
corresponds to Jt = 15. Though the field is strong; 1.8 ■ 
lOV/m (C = 6.46 ■ 10"^ ; D = 1.07), D value is similar to 
the one of Fig. [5] and the curves are correspondingly similar. 

of the pulse. Physically, the dimensionless kick strength, 
P equals to the typical amount of angular momentum (in 
the units of h) supplied by the pulse to the molecule. For 
example, in the case of KCl molecules, P = 25 corre- 
sponds to the excitation by 2ps (FWHM) laser pulses 
with the maximal intensity of 5.8 • IQ^'^W/crn?. The 
distribution functions for kicked molecules are shown in 
Fig. [51 The kick parallel to the deflecting field increases 
the rotational energy of the molecules (i.e. makes them 
untrapped), while keeping unchanged the value of the 
azimuthal momentum. As was explained in the previ- 
ous section, untrapped molecules with relatively low az- 
imuthal momentum contribute to the negative shift of 
the peak in the Ai distribution function (Fig. (8^)- Since 
most of the molecules became untrapped, the rainbow 
around 0.5 in the distribution of A2 becomes the dom- 
inant rainbow. For the numerical example under con- 
sideration, the distribution of the deflection angles has 
the same shape as the distribution for Ai (with a proper 
scaling). As follows from Fig. [8^ , a prealigning laser 
pulse applied parallel to the direction of the deflecting 
field leads to a dramatic narrowing in the distribution of 
the scattering angles, and increases the brightness of the 
molecular beam deflected by a static electric field. 

In the case of an aligning pulse in the x direction (per- 
pendicular to the deflecting field), both Pg[Q) and /^(O) 
are replaced by: 

Pq{Q) P^(0) + P;cos2(/)(0)sin(26l(0)) 

P^(0) ^ P;(O)-P;sin2(0(O))sin(2(/.(O)) (24) 

The distribution functions for Ai^ in this case are 
shown in Fig. [9l Such a pulse forces the molecules to 
rotate preferentially in the planes containing the x axis. 
In a previous work [2^ . we showed that the distribution 




FIG. 8: The statistical distributions of (a) Ai and (b) A2 for 
prealigned KCl molecules (P = 25). The prealigning pulse 
is parallel to the deflecting field (along z-axis). The deflection 
field is strong: 1.8 ■ lOV/m, and Jt = 5 (C = 5.81 • 10"" ; 
D = 9.62). 




FIG. 9: The distributions of (a) ^1 and (b) A2 for prealigned 
KCl molecules. The conditions are similar to those of Fig. 
|8l except the polarization of the prealignment pulse, which is 
perpendicular to the deflecting field. 

of A2 has two rainbows at and 0.5 in the absence of 
the static field. In the present case, as can be seen from 
Fig. inb, the distribution of A2 still preserves the rain- 
bow at 0.5, but the rainbow at is smeared due to the 
effect of the strong deflecting field {D = 9.62). Neverthe- 
less, we still may observe a considerable concentration of 
molecules at low A2 values. The peak at A2 ~ 0.28 is 
still present due to the reasons explained in the previous 
section. The two rainbows are due to molecules that ro- 
tate with small azimuthal momentum. The distribution 
of Ai (Fig. [5^) has now a strong peak at low positive Ai 
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values, since a great portion of the molecules is freely ro- 
tating with large azimuthal momentum provided by the 
laser kick. 

Finally, before proceeding to the quantum treatment 
of the same problem, we refer the reader to Sec. IVIl in 
which an alternative approach to the classical calculation 
of Ai^2 is given by means of the formalism of adiabatic 
invariants. 



IV. QUANTUM TREATMENT 

For a more quantitative treatment, involving analysis 
of the relative role of the quantum and thermal effects on 
one hand, and the strength of the pre-aligning pulses on 
the other hand, we consider quantum-mechanically the 
deflection of a linear molecule described by the Hamilto- 
nian: 



P 1 

H= iiFcose~-F\(c 

21 ^ 2 



I - a±_) cos^ 6 + aj_] , (25) 



where J is the angular momentum operator. 

Without the electric field F, the eigenfunctions of the 
free-space molecule are given by the free-rotor eigenfunc- 
tions I J, m) . Before the molecules enter the deflecting 
field, we prealign them by a short femtosecond laser 
pulse. Such a pulse creates a rotational wave packet 
of the \J,m) states. After the prealigning laser pulse 
is over, the molecules enter adiabatically the region of 
the static field. Then, each | J, m) state (within the 
wave packet) transforms into the corresponding eigen- 
state I J, to), where J is the quantum number associated 
adiabatically with the quantum number J. The relation 
between |J,to) and the free-rotor eigenfunctions may be 
described by: 



I J, to) = ^ PiJJ.m) 

J=\m\ 



The force, J- acting on the molecule is given by: 



(26) 



(27) 



The first derivative is obtained by the means of the 
Hellman-Feyman theorem, that is being in an eigenstate 
|J,to), 



= (J,to| — |J,to). 



(28) 



From Eqs. (|25]), (|27l) and ^ jt is clear that the de- 
flection angle of a molecule in a | J, to) state is given by 
Eq. (|3]), in which Ai^2 are replaced by: 



J.m 



A. 



(J, to| cos 9\J, to) 
(J, to| cos^ 9\J, to) 



In the quantum case, the continuous distribution of the 
angles 7 is replaced by a set of discrete lines, each of them 
weighted by the thermal population of the state | J, to) . 



Fig. [To] shows the distributions of Af'™ in the thermal 
case (i.e. without prealignment). These results can be 
compared to their classical analogs in Fig. [71 where the 
same structure is seen. 
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FIG. 10: Quantum distributions of (a) ^1 and (b) A2 for a 
thermal beam of KCl molecules in a strong field (1.8-loV/m) 
and Jt = 15 (C = 6.46 ■ lO"'^ ; D = 1.07). These graphs are 
similar to the classical graphs of Fig. [T] 

If the molecules are subject to a strong femtosecond 
pre-aligning pulse parallel to the deflecting field, the cor- 
responding interaction potential is given, as in the pre- 
vious section, by Eq. (|T]), in which F is replaced by the 
envelope e of the femtosecond pulse (including the 1/2 
factor, as was explained in the previous section). If the 
pulse is short compared to the typical periods of molecu- 
lar rotation, it may be again considered as a delta-pulse. 
In the impulsive approximation, one obtains the follow- 
ing relation between the angular wavefunction before and 
after the pulse applied at i = [3l| : 



*(i = 0+) =exp(jPcos2 6')*(< = 0"), 



(29) 



(30) 



where the kick strength, P is given by Eq. (|23l) . For the 
vertical polarization of the laser field, to is a conserved 
quantum number. This allows us to consider the exci- 
tation of the states with different initial to values sepa- 
rately. In order to find '^'{t — 0+) for any initial state, we 
introduce an artificial parameter ^ that will be assigned 
the value ^ = 1 at the end of the calculations, and define 

*5 = exp [{iP cos^ e)(\ ^{t = 0-) = ^cj(^)|J,to). 

(31) 

By differentiating both sides of Eq. (|?T|) with respect to 
^, we obtain the following set of differential equations for 
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the coefRcients cj: 

cj' =iP^cj{J',m\cos^6\J,m), (32) 
J 

where c — dc/d^. The matrix elements in Ea. p2p can 
be found using recurrence relations for the spherical har- 
monics [l^. Since ^^^o = ^(^ = 0") and = ^{t = 
0+) (see Eq.(|3T|)). we solve numerically this set of equa- 
tions from ^ = to ^ = 1, and find '^{t — 0^). In order to 
consider the effect of the field-free alignment at thermal 
conditions, we repeated this procedure for every initial 
I Jo, mo) state. To find the modified population of the 
|J, m) states, the corresponding contributions from dif- 
ferent initial states were summed together weighted with 
the Boltzmann's statistical factors: 

f(A \ - V- eM-E'''/kBT) 

J l-^J,mo ) ~ / , ^ 

y ^rot 

X \cj\'Sa,,^,,Aj^^^, (33) 

where cj are the coefRcients (from Eq. [32)) of the wave 
packet that was excited from the initial state | Jo, wo); S 
is the Kronecker delta symbol, and Qrot is the rotational 
partition function. The distribution in the case of parallel 
pre-alignment is given in Fig. 1111 The results are quite 
similar to the classical results from Fig. [5] 




FIG. 11: Quantum distributions for prealigned KCl 
molecules. Here the prealignment pulse (P = 25) was par- 
allel to the deflecting field (strong field, 1.8 • lO^V^/m; C — 
5.81 • 10~* ; D = 9.62). Temperature corresponds to Jt = 5. 
These graphs are similar to the classical graphs of Fig. [S] 

In the case of an aligning pulse in the x direction, the 
operator in Eg. (150)) becomes: 

^^f{t = 0+) = exp(iPcos2 0sin2 6')*(t = 0~), (34) 

and a procedure similar to the described above is used to 
find the deflection distribution (one should pay attention 



that m is no longer a conserved number during the op- 
eration of the pulse in the x direction). The distribution 
for the case of perpendicular prealignment is given in Fig. 
[T^ These results are similar to the classical predictions 
from Fig. [9| 




FIG. 12; Quantum distributions of KC'l molecules prealigned 
by the means of a pulse (P = 25) perpendicular to the de- 
flecting field (strong field: 1.8 ■ lOW/m; C = 5.81 ■ 10"" ; 
D = 9.62). Temperature corresponds to Jt = 5. These 
graphs are similar to the classical graphs of Fig. [9] 



V. DISCUSSION AND CONCLUSIONS 

In this work we considered molecular deflection by 
weak and strong inhomogeneous static electric fields. As 
the deflecting fleld is increased, it modifies the time- 
averaged alignment/orientation of the molecules. This 
affects the dipole force, and we have studied both clas- 
sically and quantum mechanically the resulting deflec- 
tion process. We found that laser induced field-free 
pre-alignment provides an effective tool for controlling 
molecular deflection. Depending on the polarization of 
prealignment pulse, different control actions may be ex- 
erted. In particaular, we predict a dramatic increase in 
the brightness of the scattered molecular beam, when the 
prealignment pulse is parallel to the direction of the de- 
flecting fleld. Though we discussed (for simplicity) linear 
molecules in this work, a similar control mechanism may 
be considered for polyatomic molecules with more com- 
plicated geometry. Being in free space such molecules 
rotate about their own axis as well, which leads to A2 
distribution different from that of Eq. ([7]) , and Ai is not 
necessarily 0. For such molecules the prealignment will 
play a signiflcant role in reducing their dipole interaction 
with the static fleld. Molecular deflection by inhomoge- 
neous static electric flelds may be used for the separation 
of molecular mixtures. Narrowing the distribution of the 
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scattering angles may substantially increase the efficiency 
of separation of multi-component beams, especially when 
the pre-alignment is applied selectively to certain molec- 
ular species, such as isotopes fsS^, or nuclear spin iso- 
mers 34, 35] . Controlling molecular deflection by means 
of laser-induced prealignment may be implemented also 
for magnetic molecules moving in the static inhomoge- 
neous magnetic fields, and this phenomenon is a subject 
of the currently ongoing research. Controlling the dipole 
interaction by laser-induced pre-alignment may find ap- 
plications in molecular deceleration methods using time 
and spatially varying electric and magnetic fields [36]. 

We acknowledge the support of our study by a grant 
from the Israel Science Foundation. This research is 
made possible in part by the historic generosity of the 
Harold Perlman Family. lA is an incumbent of the Pa- 
tricia Elman Bildner Professorial Chair. 



VI. APPENDIX: CALCULATION OF Ai,2 BY 
MEANS OF THE THEORY OF ADIABATIC 
INVARIANTS 

The energy of a molecule participating in the deflection 
process is: 

H ^^I (i)^ + <j>^sm^0^ - nFcos9, (35) 

where we neglected the effect of polarizability (which is 
small for the fields and molecules considered in the main 
body of the paper). The conjugate momenta and Pg 
are given by Eqs. ([T4|) and ([T5|) . respectively, and is a 
constant of motion. It is convenient to change variables 
and to define new constants [13, : 

u=— cos 9, (36) 



where g{u) is the rhs of the Eg. ([39]) . If a = 0, the 
polynomial g{u) has two roots, ui,U2 such that —1 < 
ui < U2 < 1. When a ^ 0, g{u) has three roots (de- 
noted by Mi,U2,U3). Let us analyze the behavior of this 
polynomial. If u 3> 1, then g{u) « av? > 0. If we 

substitute u = 1, then g{l) = - (^^"j < 0. That 

means that the three real roots of g(u) are ordered as 
— 1 < ui < M2 < 1 < M3. The polynomial g{u) is positive 
between ui and U2, and each root corresponds to a real 
angle 9. This allows one to determine the half-period of 
the motion in the static field F by integrating Eq. pO)) 
from 1*1 to U2. 

Since the potential is time-dependent (at least dur- 
ing the rising time), the energy of the system is not a 
constant of motion. However, the defiection potential 
is adiabatic with respect to the rotational motion and, 
therefore, we can use adiabatic invariants to determine 
the energy of the system [37l - [39| . The adiabatic invariant 
related to the coordinate 9 is: 

Ig= Ped9. (41) 
From EQS. (fT5|) . and it is easy to derive: 

where k is a constant. The energy H of the molecule 
as a function of the energy Hq without electric field is 
obtained numerically by solving the following equation: 



/? = jH, (37) 



le = II (43) 



where Ig is calculated for a = 0, that is in the absence 
of the external field. 

2u,F 

a = — — . (38) Once the energy of the system H and the polynomial 



g{u) are found, Ai^2 is given by: 



I ' 

It is easy to show that u obeys the following equation: 



which gives: 



/77/ 

dt=^=, (40) Eq.dH 



This is equivalent to calculating Ai,2 with the help of 
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